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Abstract 

We propose a novel dipole trapping scheme using 
spatial light modulators (SLM) for the manipu- 
lation of individual atoms. The scheme uses a 
high numerical aperture microscope to map the 
intensity distribution of a SLM onto a cloud of 
cold atoms. The regions of high intensity act 
as optical dipole force traps. With a SLM fast 
enough to modify the trapping potential in real 
time, this technique is well suited for the con- 
trolled addressing and manipulation of arbitrar- 
ily selected atoms. 

1 Motivation 

Many proposals for large-scale quantum comput- 
ing [JJ or quantum simulation [2] rely on the abil- 
ity to detcrministically manipulate, address and 
couple individual components of a quantum net- 
work. These challenges can be adressed by the 
use of isolated quantum systems, such as trapped 
ions in quantum computing and entanglement ex- 
periments H2E], or dipole-trapped single neutral 
atoms [SHI]- In an effective one-dimensional ar- 
rangement it has already been shown that optical 
addressing and preparation of individual dipole- 
trapped atoms is possible and that those atoms 
can also be used as a phase preserving quan- 
tum register [5]. To control a larger number 
of atoms in a two or three-dimensional quan- 
tum network, more elaborate techniques are re- 
quired. Therefore several methods for handling 



individual atoms have been developed. The 
most promising techniques are the trapping of 
ion strings [21I31H], magnetically trapped neu- 
tral atoms above atom chips [TTTMl2j . individual 
atoms in steep optical dipole-force traps [6 , opti- 
cal lattices [7] , and dipole-trap arrays created ei- 
ther by a matrix of lenses [12] , by holograms Q3] , 
or by a combination of these techniques |15) . 
In contrast to these attempts, we propose to 
use a spatial light modulator (SLM) [TB1[T7] to 
form the desired trapping potential. This con- 
cept has much in common with the optical tweez- 
ers technique used in biology [TH], where light is 
used to manipulate small objects, like viruses, 
bacteria, or organic samples attached to micro- 
spheres. Here we propose to apply this technique 
directly to individual atoms. This combination 
of the SLM with optical tweezers offers atom- 
chip like flexibility of individual-atom manipula- 
tion as well as magnetic-field free dipole trapping. 
Furthermore, with SLM refresh rates far above 
the loss rate of trapped atoms, such an arrange- 
ment might allow for the real-time manipulation 
of atoms. 



2 Outline 

At the heart of this novel scheme is the use of a 
SLM to generate arbitrary patterns of light in 
the object plane of a microscope and to trap 
atoms therein. Fig.[T] illustrates this approach 
and shows how one can use the same microscope 



to observe the atoms. Several types of SLM are 
available which all have the potential to achieve 
this goal. One widely used technology is the liq- 
uid crystal device (LCD) [T3JQJ]. It has the ad- 
vantage of creating a range of grey scales, but has 
the drawback of having a low refresh rate, usu- 
ally not exceeding 100 Hz. This is in general too 
slow for the dynamic control of trapped neutral 
atoms. Other commercially available SLMs are 
digital mirror devices (DMD) [20] . They consist 
of an array of individual flat mirrors which can 
be independently switched between two tilt an- 
gles to generate arbitrary intensity distributions. 
With a full-frame refresh rate of up to 50 kHz 
and a resolution of typically 1024 by 768 indi- 
vidual mirrors, these devices seem ideal for the 
real-time manipulation of trapped atoms. 

The desired intensity distribution can be 
formed in two ways. One possibility is to place 
the SLM in the back focal plane of a lens sys- 
tem, such that its Fourier transform is in the 
front focal plane. Thus the SLM acts as an ef- 
fective hologram for the resulting intensity distri- 
bution [T3] . This method has the drawback that 
any change in the trapping pattern requires one 
to recalculate the entire hologram and thus re- 
strains the experimenter to work with predefined 
sequences. Alternatively, the tweezers can be 
generated by imaging the surface of an amplitude- 
modulating SLM through a microscope directly 
onto the atoms. The tweezers are then formed 
in the object plane of the microscope, see Fig. [2] 
In the following we concentrate on this latter 
method, as the direct mapping of intensities pro- 
vides an intrinsically higher speed and flexibility 
in the manipulation of trapped atoms. This sys- 
tem can be used to expose the atoms to almost 
any arbitrarily shaped and time-varying poten- 
tial landscape, which offers a plethora of applica- 
tions. Out of these we focus on the implementa- 
tion of an array of very tiny atom traps, eventu- 
ally capable of storing single atoms. 

Any tight spatial confinement of atoms in 
dipole-force traps relies on small foci of the trap- 
ping light. Therefore, an optical system of high 
numerical aperture is necessary. This can be 
done with either a sophisticated microscope ob- 
jective [21] or a single aspheric lens [22]. Typi- 
cally, these are diffraction limited and have a nu- 
merical aperture of 0.5, leading to an optical reso- 




Figure 1: Scheme of the proposed optical tweez- 
ers. A red detuned laser beam illuminates the 
mirror array of the SLM. A microscope repro- 
duces the intensity pattern of the SLM at the sur- 
face of a mirror placed in the tweezer plane (ob- 
ject plane of the microscope) . The reflected light 
interferes with the incoming beam and forms a 
standing wave. Its antinodes act as dipole traps, 
which can hold individual atoms. Since the SLM 
is reconfigurable the setup allows arbitrary repo- 
sitioning of atoms. A very sensitive CCD camera, 
detecting atomic fluorescence via a dichroic mir- 
ror, can be used to image the trapped atoms and 
monitor the tweezers in operation. The insert 
(Courtesy of Texas Instruments) shows a block 
of 3 by 3 mirrors of the SLM. The mirrors can be 
switched between two tilt angles, so that the light 
is directed either into the microscope or deflected 
to a beam stop. 



lution approximately equal to one wavelength, A. 
This also determines the lateral confinement of 
the atoms. Along the optical axis, the atoms are 
trapped by the longitudinal beam profile. For a 
single beam focussed to A, atoms are axially con- 
fined to a Rayleigh length of ir-\. To obtain a bet- 
ter confinement along the optical axis, a stand- 
ing wave could be used. This can be achieved by 
placing a mirror at or close to the plane where the 
tweezers are formed. The atoms trapped in the 
anti-nodes of the standing wave are then axially 
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Figure 2: The top picture shows a pattern on 
the SLM. The pattern is a repetition of a 4 by 
4 mirror block where the corner pixels are miss- 
ing. Black marks the mirrors which are switched 
on. Each mirror has a dimension of 14 /im by 
14 /im. The bottom picture shows the resulting 
image, where black is the highest intensity, when 
an imaging system with a numerical aperture of 
NA = 0.5 and a magnification of 38 is used. 

The underlying physical mechanism used to 
trap atoms in optical tweezers is the optical 
dipole force, which atoms experience in detuned 
light fields. This force is conservative but very 



weak, resulting in a shallow potential with a typ- 
ical depth of ~ 1 mK. Only atoms with low ki- 
netic energy can be trapped, and therefore they 
need to be pre-cooled before loading them into 
the tweezers. This is usually accomplished by a 
magneto-optical trap (MOT), trapping and cool- 
ing atoms down to the Doppler cooling limit, 
To ~ 100 /zK for most alkali atoms 23 . In 
the case where a mirror is included to form a 
standing-wave dipole trap, the pre-cooling has to 
be done close to the mirror surface. This can 
easily be accomplished using a magneto-optical 
surface trap (MOST) 

3 Feasibility 

3.1 Dipole trapping of Alkali 
Atoms 

The majority of experiments cooling and trap- 
ping neutral atoms are done with alkali atoms, 
of which Rubidium is one of the most prominent 
species. Hence we choose it for the purpose of 
this feasibility study. The optical tweezers con- 
sist of tiny optical dipole-force traps, in which 
the dynamic Stark shift gives rise to a trapping 
potential in the far detuned limit [2 5) : 



Udi P (r) 



hT 2 J(r) 
85 I sat 



(1) 



where 8 is the detuning of the trapping laser light 
with respect to the atomic transition, I(r) is the 
intensity of the trapping light, I sa t is the sat- 
uration intensity for the chosen transition, and 
T is the corresponding decay rate. For red de- 
tuned light, 5 < 0, Udip is negative and hence 
attractive. The most critical parameter is the 
trap depth Uo, which corresponds to the highest 
intensity of the trapping laser light Iq. Hence 
Uq is the energy required for an atom at rest to 
escape the trap, neglecting gravity. Obviously, 
the trap depth has to be at least as large as the 
energy of the pre-cooled atoms. Therefore, the 
minimum usable trap depth is determined by the 
temperature that can be reached by the magneto- 
optical trap. For rubidium, the Doppler cooling 
limit would be To = 143 /iK. However, it is de- 
sirable to operate deeper traps since the spatial 
confinement of the atoms increases with depth. 
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For atoms whose kinetic energy is smaller than 
the potential depth by a sufficient amount, the 
trap can be treated as harmonic, and the spatial 
confinement of the atoms can be determined by 
the harmonic oscillation frequencies u r and u> z 
for radial and longitudinal motion, respectively. 
The achievable values of Uq, ui r , and oj z depend 
on the wavelength A of the trapping laser, its in- 
tensity, and the trap geometry. 

In the proposed scheme (Fig. [I]), the trap ge- 
ometry is essentially determined by the distance 
L between the focal plane of the optical system 
which generates the tweezers and the dielectric 
mirror. The intensity modulation induced by the 
interference of forward and backward travelling 
waves over this distance can be seen in Fig. [3] 
In the following, we first discuss the limiting case 
L — > oo (left picture of Fig. |3J), which corre- 
sponds to having no mirror at all, and then an- 
alyze the longitudinal intensity modulation for 
finite L. 

When the dielectric mirror is not present (L — > 
oo), the trapping potential has a single minimum 
at the beam waist, in the focal plane of the lens 
system. The maximum achievable intensity at 
that point, 

I a = I d -M 2 , (2) 

depends of the magnification M of the optical 
system and the damage threshold the DMD with 
respect to the incident intensity, which is about 
Id = 10 W/cm 2 . To maximize the spatial control 
of the atoms, the individual traps should be as 
small as possible. The smallest achievable trap 
size is set by the resolution of the optical system 
D m i n . According to the Raylight criteria it is: 

D min = 1.22 • XN, (3) 

where A is the wavelength of the light used, and 
N is the f-number related to the numerical aper- 
ture NA of the optical system by N « ^Wa ■ ^ an 
aspherical singlet lens is used, the NA can be as 
high as 0.5. For this value, the resolution limit 
is given by about the wavelength of the trapping 
light, which is A s» 780 nm for the D2 line of Rb. 

3.2 The role of the SLM 

It is important to bear in mind that any experi- 
ments involving dynamics require the reconfigu- 
ration of the trapping potential. One way atomic 
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Figure 3: Intensity distribution of a focused laser 
beam (with beam waist wq = 0.78 /im) in the 
neighborhood of the focal plane for different dis- 
tances L between the dielectric mirror and the 
focal plane. (Black corresponds to the highest 
intensities.) The left picture is for L — > 00, i.e. 
no mirror being present. The tweezer is a pure 
forward traveling wave with no longitudinal in- 
tensity modulation. The middle picture shows a 
beam which is reflected in its focal plane (L = 0). 
The tweezer has a perfect standing wave pattern, 
but atoms cannot be trapped in the beam waist. 
The right picture shows the intensity of a retro- 
reflected beam whose focus lies L = 8/im above 
the mirror surface. Here the intensity modula- 
tion is notably less than in the previous case. 



transport could be realized is by displacing a pat- 
tern of traps. The switching of mirrors would 
lead to abrupt changes in the trapping poten- 
tial if the contribution of light for a point in 
the image plane of the DMD is dominated by 
one mirror only. If, however, the optical system 
does not resolve individual mirrors, the reconfig- 
uration of the traps can be smoother, since the 
light intensity for every point will have contri- 
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butions from a number of mirrors. Effectively, 
this leads to grey scales in the trapping pattern. 
To take advantage of this, we choose the mag- 
nification M such that blocks of 2 x 2 mirrors 
will just about be resolved. The micro-mirrors 
of a commercially available DMD have a size of 
14 fxm x 14/im each, so in our case, a magnifica- 
tion of M — 38 is required for a diffraction limit 
of D min m 780 nm (see Fig. 0). According to Eq. 
©, if the DMD is illuminated with light of an 
intensity at its damage threshold, an intensity Iq 
of 14.4kW/cm 2 can be obtained at the center of 
a tweezer. The saturation intensity for the D2 
line of Rubidium is I sat — 2.5mW/cm 2 . There- 
fore the maximum intensity is 5.8 • 10 6 times the 
saturation intensity. Taking the above discussion 
into account and requiring a trap depth of 1 mK 
leads to a detuning S < 2.1 • 10 5 r, according to 
Eq. fl]). This corresponds to a wavelength de- 
tuning of about 2.6 nm to the red with respect 
to the D2 line. However, since Rubidium has a 
fine structure splitting of 15 nm, this estimation 
of the required detuning is certainly too crude. 
The hyperfinc transitions within the Di and D2 
lines of Rb need to be taken into account. Two 
possible wavelength regimes present themselves. 
The first is slightly red detuned from the D2 line 
at 780 nm, and the second is red detuned from 
the Di line at 795 nm (see the top graph of Fig. 
Q|. More specifically, for a trapping laser with 
an intensity of Jo = 14.4 kW/ cm , for instance, a 
trap depth of 1 mK can be achieved for the wave- 
lengths of A 2 = 782.85 nm, and Ai = 796.90 nm. 

Working with light so close to the Di and 
D2 resonnances gives rise to atomic heating due 
to photon scattering from the trapping beam. 
The heating is proportional to the scattering rate 
Rscat, which, in the far detuned limit |25j . is 
given by: 

r 3 j 

Rscat ~ 7Tf2 7 ' (4) 
8iH I S at 

Here we have to take the scattering from both 
the DI and the D2 line into account, as shown 
in the bottom graph of Fig. 0} For both Ai and 
A2, Rb scatters about 600 photons per second. 
The scattering rate in both detunings is about 
the same, since in both cases, the detuning to 
the closest transition is small when compared to 
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Figure 4: The top graph shows the potential 
depth of the trap with respect to the wavelength 
of the trapping light, for linearly polarized light 
with an intensity of Iq = 14.4kW/cm 2 . The 
trapping regions are hatched. The bottom graph 
shows the corresponding scattering rate. 



detuning from the other transition. If the trap- 
ping laser were further detuned, the region to the 
red from the Di line would be preferable, as the 
potentials of both lines add constructively there, 
whereas in the region in-between the Di and D2 
lines, the potentials add destructively. 

Given that each scattering event increases the 
kinetic energy of an atom by an amount equal 
to the recoil energy, the lifetime of a Rubidium 
atom (initially at rest) in the dipole trap would 
be 2.3 s if photon scattering was the only source 
of heating. It will be shown below that this time 
is long enough to allow for many transport and 
detection cycles. 
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3.3 Radial and axial confinement 

We now estimate the radial and longitudinal con- 
finement for a travelling-wave dipole trap {L — > 
oo). The intensity of the light field in the radial 
direction is 

where q = yjL r is the distance from the axis 
of the trap, and TV sa 1/(2 NA) is the f-number 
which is about 1 in this case. Ji{q) is a Bessel 
function of the first kind. For small q, it expands 
as 

Ji(q) = \q-^i 3 + o(q 5 ). (6) 

Inserting this into ([5]) and keeping only terms up 
to second order we obtain the harmonic approx- 
imation to the trapping potential in the radial 
direction 

where Uq is the trap depth. It is plotted in Fig. 
Has a function of the wavelength A. A particle 
with mass m oscillates with frequency 

^ = V^A^ (8) 

in this potential. For — Uq = 1 mK ■ ks, this 
leads to a radial trapping frequency of co r = 
27r • 162 kHz. If we approximate this as a Gaus- 
sian distribution, the beam waist would be wq = 
V2XN/tt = 0.61 fim, which is on the same or- 
der of magnitude as the resolution limit deter- 
mined with ([3]). The oscillation frequency along 
the axial direction uj z can be calculated using the 
Rayleigh length zq = ix Wq/X « 1.5 /zm: 

- - « 9 > 

For — Uq = 1 mK • ks , this yields w z = 2tt ■ 46 kHz. 
The axial trapping frequency lu z is significantly 
lower than the radial trapping frequency u) r , 
which is generally the case for running wave op- 
tical dipole- force traps. 

We now consider the dielectric mirror to be in 
the object plane of the microscope L = 0. The 
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Figure 5: Potential depth U(z) along the axis, 
with the focal point 8 fim away from the mirror 
surface. The slowly varying average of the po- 
tential (gray) is a feature of focusing the light, 
whereas the oscillation (black) is caused by the 
partial standing wave. The modulation has the 
highest visibility close to the mirror-surface, be- 
comes less pronounced in the focal spot, and then 
increases again further away from the focal point. 



backward travelling wave interferes with the for- 
ward travelling wave and gives rise to a longitudi- 
nal interference pattern, creating an ideal stand- 
ing wave (see the central picture of Fig. 3). The 
potential has many local minima in which atoms 
can be trapped. The intensity modulation along 
the z-axis gives rise to the following longitudinal 
potential: 

U(z) = A(z)sin 2 . (10) 

The envelope of the modulation, A{z) = 
4Uq/ (1 + z 2 /zq), reaches its maximum at the mir- 
ror surface. Here Uq is the potential depth only 
due to to a single focused beam. Because the 
contrast of the modulation is one, A(z) is also 
the trap depth of a local minimum at height z 
above the mirror. The largest achievable trap 
depth is 4[/o, four times higher than without 
mirror. Therefore the maximum radial trapping 
frequency in the standing wave configuration is 
wp" = 2ui r = 2ir ■ 324 kHz and is achieved close 
to the mirror. The longitudinal oscillation fre- 
quency inside a local trap at height z is given 
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by 



-8A(z)ir 2 



i\ 2 



ill) 



Close to the mirror, the highest value wf* = 2tt ■ 
1100 kHz is achieved. 

When the mirror is not exactly in the image 
plane of the DMD (finite distance L), the for- 
ward and backward propagating waves do not 
have the same amplitude. Therefore the longitu- 
dinal intensity modulation, whose peak-to-peak 
amplitude is 

4t/ 



A(z) 



V'l + zVz Vl + (2i-z) 2 Ao 2 ' 



(12) 



rapidly fades away as the distance between the 
mirror and the focal plane increases. Fig. 5 
shows the longitudinal potential U(z) for L = 8 
/im. The deepest local minima are found at the 
focus of the beam. An atom trapped in one of 
these local minima sees a trap depth of about 
1.37 Uo in the radial direction but only 0.3 Uq 
in the longitudinal direction. Whatever the dis- 
tance between the focal plane and the mirror, 
the trap depth in the radial direction will al- 
ways be larger than in the travelling wave case 
(L — > oo) but smaller than in the standing wave 
case (L = 0): 162 kHz < w r /(27r) < 324 kHz. 
Along the z-axis, the trapped atoms are con- 
fined cither to the slowly varying envelope or to 
the local minimum, depending on their initial 
kinetic energy. Only the coldest atoms can be 
trapped in the local minima. If they are, their 
confinement shall be strongly improved. In that 
case the longitudinal frequency can be estimated 
from Eq. QTJ with A(z) given by Eq. Q21 For 
the example of Fig. 5, one finds that the os- 
cillation frequency at the focus of the beam is 
lu z = 2tt ■ 300 kHz. Depending on L, the values of 
lo z are intermediate between the travelling wave 
case (L — > oo) and standing wave (L = 0) case: 
46 kHz < w z /(2tt) < 1100 kHz. 

4 Ways to single atoms 

In addition to controlling the size and depth of 
the traps, it is also important to control and 
measure the number of atoms. More specifi- 
cally for the goal of conducting single atom ex- 
periments an efficient scheme to prepare single 



trapped atoms is essential. When loading atoms 
into moderately sized traps, the number of atoms 
therein is generally Poisson distributed. In this 
case, we would need to measure the atom num- 
ber and then post-select traps with just one atom. 
Fortunately, for tightly focused dipole traps, the 
collisional blockade mechanism [27] should give 
rise to a significant departure from the Poisso- 
nian statistics and favor the loading of only one 
atom per trap. 

The effect occurs when atoms are loaded at a 
rate R from a magneto optical trap into tightly 
confined dipole traps. It relies on light-assisted 
collisions [25] between atoms in the presence of 
the MOT cooling beam. Provided the latter is 
strong enough to saturate the transition, two- 
body losses occur at a rate {3'N(N — 1). N is the 
atom number and /3' is the rate constant, which 
is inversely proportional to the trapping volume. 
There is also a single body loss process 7 TV, which 
is mostly due to background collisions, so the full 
rate-equation [29 for the entire loading process 
reads: 



dN 
~dt 



R-jN - P'N(N - 1). 



(13) 



Collisional blockade will occur when the two- 
body loss rate dominates over the loading rate 
of the dipole trap. Hence, it will work very ef- 
ficiently when f3' is large, i.e. if the atoms are 
confined to a tiny volume. The volume occu- 
pied by the atoms depends on the temperature 
and the trapping frequency. A comparison of the 
frequencies achievable with the proposed set up, 
and the trapping frequency io t = 2ir ■ 200 kHz 
of a previous experiment demonstrating the col- 
lisional blockade [27] leads us to the conclusion 
that this mechanism could here be highly efficient 
as well. 

5 Single atom detection and 
transport 

To investigate the loading processes, the atom 
numbers in the individual traps have to be de- 
termined. For well-isolated traps, laser induced 
resonance fluorescence could be used for atom 
counting [26. . To collect the fluorescence, the lens 
system which creates the optical tweezers can be 
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used in the backward direction as a microscope. 
This fluorescence is detected by a highly sensitive 
camera. 

To discuss the feasibility of single atom de- 
tection, the flux of photons impinging on the 
camera has to be determined. The geometri- 
cal collection efficiency 77 depends on the numer- 
ical aperture of the lens system. To evaluate 
77, we assume a uniform photon emission of the 
atoms and calculate the ratio of photons pass- 
ing through the system. The ratio of the sur- 
face area of a spherical cap with that of a sphere, 
yields i] = 2n (l — cos(f )) /(47r), where a is the 
opening angle. For a numerical aperture of 0.5, 
a = 60°, we expect a collection efficiency of 
rj = 0.067. 

If the radiative transition of the atom is com- 
pletely saturated, the scattering rate on reso- 
nance tends to R sca t — |r, hence R sca t = 
19.6 /is -1 for Rubidium. Assuming a further loss 
of 50% through the lens system, a camera should 
collect photons from a single atom with a rate 
of Ream = 653 ms -1 . In other words, a cam- 
era operating with an exposure time of 100 /is 
could count as many as 65 photons. For mod- 
ern electron multiplying charged coupled device 
(EM-CCD) cameras, this is well above the signal 
to noise threshold, so it is possible to use such 
devices to observe single atoms. 

The high frame rate of the DMD should allow 
for the reconfiguration of the potential landscape, 
and hence for the transport of confined atoms. A 
naive method would be to transport atoms adia- 
batically by changing the potential slowly enough 
so that the atoms are following without being 
heated. Since the DMD has only a finite res- 
olution, this must be done using discrete steps. 
With the aforementioned trap consisting of a four 
by four block, as shown in Fig. ® the pattern 
would have to be shifted by at most one pixel 
per step. The longest distance the traps could 
be moved is 1000 pixels, that is from one edge 
of the DMD to the other. With a refresh rate of 
50 kHz, this transport would take 20 ms. This is 
much shorter than the anticipated lifetime of the 
trap, which should be of at least several 100 ms. 
However a drawback of this method is that dis- 
crete switching of mirror elements leads to para- 
metric heating of the atoms. To circumvent this, 
the atom-transport could also be realized ballis- 



tically. The initially trapped atoms would be ac- 
celerated through a sudden change in trapping 
potential, guided through a channel and eventu- 
ally decelerated and recaptured. This scheme has 
the advantage that it needs only a few changes in 
the potential, making it less susceptible to heat- 
ing, and yet faster at the same time. 

6 Conclusion and outlook 

Our study shows that the handling and manipu- 
lation of individual dipole-trapped atoms should 
be possible with an optical-tweezers type setup 
with the technologies at our disposal today. Us- 
ing this new scheme, real-time arrangement of 
dipole traps is possible, and two-dimensional 
transport of single atoms seems feasible. We ex- 
pect this technique to boost the field of quan- 
tum computing, as it will allow to handle sin- 
gle entities within a large scalable array. For in- 
stance pairwise entanglement through either cav- 
ity QED [30], or controlled collisions [3T] could 
be used for generating large cluster states which 
are the essential resource for one way quantum 
computing (32] . 
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